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We study the dynamics of entanglement transfer in a system composed of two initially correlated 
three-level atoms, each located in a cavity interacting with its own reservoir. Instead of tracing out 
reservoir modes to describe the dynamics using the master equation approach, we consider explicitly 
the dynamics of the reservoirs. In this situation, we show that the entanglement is completely 
transferred from atoms to reservoirs. Although the cavities mediate this entanglement transfer, we 
show that under certain conditions, no entanglement is found in cavities throughout the dynamics. 
Considering the entanglement dynamics of interacting and non-interacting bipartite subsystems, 
we found time windows where the entanglement can only flow through interacting subsystems, 
depending on the system parameters. 
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I. INTRODUCTION 

Entanglement has emerged as a central physical re- 
source for Quantum Information Theory [lj . For process- 
ing quantum information, physical architectures are ex- 
pected to be composed of multipartite quantum systems. 
Understanding how quantum entanglement is transferred 
between the parties has motivated several contributions 
in recent years [2|-|6J]. One issue is entanglement flow 
through individual parties and the whole multipartite 
system |2|. Another interesting problem is entanglement 
transfer between qubits and its relation to energy [3| . In 
addition, we have the study of the entanglement transfer 
between non-interacting qubits, leading to conservation 
rules for entanglement depending on how qubits are ini- 
tially correlated |4J. Entanglement transfer from atoms 
to cavity modes leading to entanglement revivals has 
been studied in [5[ . Entanglement transference from two 
cavities to their corresponding reservoirs allowed the de- 
scription of Entanglement Sudden Death as opposed to 
Entanglement Sudden Birth jgj. 

An interesting problem in this context, is entangle- 
ment flow between parties in a multipartite system, in- 
cluding their dissipative mechanisms. In this work, we 
study the entanglement dynamics of two initially corre- 
lated atoms placed in two non-interacting leaky cavities, 
each connected to its own reservoir. To achieve this, we 
developed a hybrid analytical approach for finding the 
quantum dynamics of the atom-cavity-reservoir system. 
Unlike the master equation approach, our method allows 
us to include the reservoir dynamics, thus preventing in- 
formation loss due to trace operations. We study the 
evolution of entanglement in different non-interacting bi- 
partite subsystems, such as atom-atom, cavity-cavity and 
reservoir-reservoir. We show that the entanglement ini- 
tially contained in the atomic subsystem is completely 
transferred into the reservoir-reservoir subsystem. Al- 
though cavities are the bridge connecting atoms to reser- 
voirs, we show that they may not be entangled through- 



out the dynamics. Moreover, depending on the initial 
state, quantum dynamics may led to a situation where 
no entanglement exists in any of these three subsystems. 
In this case, we extend the study to other interacting and 
non-interacting subsystems. 

This paper is organized as follows: in Section II we 
develop a hybrid analytical method to find the quantum 
dynamics. In Sec III, we study the dynamics of entangle- 
ment transfer. In Section I V we present our concluding 
remarks. 



II. ATOM-CAVITY-RESERVOIR DYNAMICS 

Our model considers two independent subsystems, 
each formed by a three-level atom inside a leaky QED 
cavity. Each atom interacts with a single mode of fre- 
quency u of the quantized electromagnetic field and a 
classical field with frequency v in a Raman configuration, 
as shown in Fig. [TJ The quantum mode couples level \g) 
and \c) while classical field couples levels |e) and |c). As- 
suming no direct coupling between cavities, the dynam- 
ics of each atom-cavity-reservoir system can be studied 
individually. Neglecting effects of spontaneous emission 
from leves \c) and e), the Hamiltonian describing this 
system can be conveniently written, in the rotating wave 
approximation, as: 
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H = hA\c)(c\ - hS\e){e\ - h^(w - uj k )b\h k 

k=i 

+m(\c){e\ + \e)(c\) + hg(a\c){g\ +a t | ff )(c|) 

N 

+h^{g k abl + g* k b k tf). (1) 
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Here, a(cfi) annihilates (creates) a photon with fre- 
quency ijj in the cavity mode, operator b k {b k ) annihilates 
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FIG. 1. Scheme of a three-level atom inside a cavity coupled 
to a reservoir. 



(creates) a photon with frequency W& in fc-th mode of 
the reservoir. We have defined A = uj cg — w, with uj cg 
is the frequency difference between levels \c) and \g), de- 
tuning 6 = uj ce — v — A. The coupling strength between 
the classical field and the atom is f2 and we have ex- 
plicitly written the system-bath interaction with a linear 
coupling with coupling constant <?&. In the usual master 
equation approach, considering the Markov approxima- 
tion and an infinite number of bath oscillators, we can 
describe the dynamics of the atom-cavity system by 



p=-i[H',p] + -(2o)pa 



a* ap 



pa* a) 



(2) 



where H' corresponds to Hamiltonian ([T} without terms 
involving the reservoir. 

In such an approach, only atom-cavity dynamics is de- 
scribed, while the reservoirs are traced out. This results 
in loss of quantum correlations between reservoirs and 
other parties. 

In what follows, we will develop a method to find the 
dynamics, including explicitly the reservoir's degrees of 
freedom. In this manner we preserve all the information 
about quantum correlations in the system 

To address this problem, we develop a hybrid analyti- 
cal approach to find the dynamics. First, we will follow 
the method described for obtaining the quantum dynam- 
ics and the entanglement properties of inhomogeneously- 
coupled systems [3, Q . This method consists of inspect- 
ing the Hilbert space occupied for the quantum system 
throughout the evolution and implementing a trunca- 
tion criteria based on a probabilistic argument. Then, 
to find explicitly analytical expressions for the quantum 
dynamics, we follow the well-known quantum trajectory 
method [nj. 

Let us consider the case of a single initial excitation 
contained in the atomic subsystem. That is, 



where, |0) c corresponds to the vacuum of the cavity and 
|0) r = fl k |0fc) denotes the vacuum of the reservoir. Now 
we have to look for the accessible Hilbert space for the 
compound system when starting from this initial condi- 
tion. We can find this by following the action of Hamil- 
tonian Hu on the initial state |V>o}- It is not difficult to 
realize that a portion of the Hilbert space connected by 
Hamiltonian (p} is given by states {|e00), |c00), glO)}. 
The one photon excitation in state \glO) is transferred 
to bath modes through Hamiltonian ((T|) as follows: 



H u \glO) = %|c00) + hN \g01o), 
where we have defined the state: 
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Lo)r = -TT- y^gfc|lfc)r- 
fe=l 



(4) 



(5) 



with No = y/Y?k=i ISfcl 2 ! an d |lfc)r corresponds to the 
state having one photon in the fc-th reservoir mode and 
zero photon in the remaining modes. The state in Eq. ([5]) 
is a collective state of the reservoir having a single exci- 
tation. The |<?01o) collective state evolves under the in- 
teraction part in Hamiltonian ([1}, back to states \glO) 
and |cOO). However, this is no longer true when consid- 
ering the action of the reservoir Hamiltonian on reservoir 
states. Taking the vacuum state of the reservoir, we have 
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X)(w-w fc )S£S fc |o) r = o. 



For |1o)d free energy term leads to 



(6) 
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^ = TT^lsk{^ - Wfc)|lfc) r . 



(7) 



Although this state is different from |lrj)r in Eq-©, it 
can be written as a superposition of this state and a state 
|li) r orthogonal to |lo)r- That is, 

|ll>r = -j^[|*l>-(lo|*l>|lo>r], (8) 

where N x = «*i|*i> - |(l |*i}| 2 ) 1/2 . 

The new generated state |li) r leads to other orthogo- 
nal states having one excitation through the action of the 
reservoir terms of the Hamiltonian [7|. In consequence, 
the accessible Hilbert space for the overall system ini- 
tially prepared in state ([3]) , can be written in a collective 
basis spanned by the set of orthogonal states 

{|eOO), |c00>, |. 9 10), | ff 01o), |flOli), | 5 0l 2 ), . . . }. (9) 

Therefore, the atom-cavity-reservoir system governed by 
Hamiltonian ((T|) evolves to 



\ipo) = \e) a <8> |0) £ 



|eOO), 



(3) 



\iln) = E t \e00) + C t \dOO) + G t \glO) 
+Ro >t \g01o) + Ri,t\g01i) + . 



(10) 



where Rjj is the time-dependent probability amplitude 
of the state |<?01j). This state can conveniently written 
defining one-excitation collective state for the reservoir 
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(11) 



(12) 



corresponds to the probability amplitude of having one 
excitation in the reservoir modes. 

In this manner, the dynamics of the atom-cavity- 
reservoir system can be described in terms of a three- 
qubit system in the basis: {|e), \g)}®{\0), |l)}<g>{|0), 1)}. 
Notice that we have assumed the high detuning regime 
(A 3> 0,,g), so the electronic level \c) a is only virtually 
populated and the evolution described in Eq. (|10[) can 
now be written as 



\il>t)=E t \eOO)+G t \glO)+Rt\g01). 



(13) 



Probability amplitudes can be obtained by numerical di- 
agonalization of Hamiltonian ([1]) where it is found that 
Ct ~ 0. It is worth to notice that this form of finding the 
quantum dynamics has close relation with the Weisskopf- 
Wigner procedure (see for example ref. [9J]) where the 
reservoir is considered throughout the dynamics. 

We are interested in finding analytical expressions for 
the temporal dependent coefficients in the wave function 
in Eq. (fT3|). Unfortunately, this is not possible by solv- 
ing the Schrodinger equation. However, if we know the 
reduced dynamics for the atom-cavity system, we can get 
information about the probability amplitude Rt- The re- 
duced atom-cavity dynamics could be obtained by solv- 
ing the corresponding master equation, or alternatively 
by using the quantum trajectory method (see for exam- 
ple [l(| and references within). By using this approach, 
we obtain the reduced density matrix for the atom-cavity 
system as: 

p(t) = \Et\ 2 \eO){eO\ + \G t \ 2 \gl){gl\ 

+E t G* t \eO)(gl\+E;G t \gl)(eO\ 

+ |i? t | 2 |<?0)( 5 0| (14) 

The probability amplitudes obtained through the quan- 
tum trajectory approach are 

E t = {cos (fit) + -£L sin (Ctt)}e~* Kt , 



G, 






(15) 



where 4ti 2 = 4<?| ff — k 2 /4 and g c ff — gfl/A. For sim- 
plicity, here we have set 6 = (g 2 — n 2 )/A. From these 
expressions we obtain the probability amplitude Rt as: 



\Rt\ = VI - \Et\ 2 - \G t 



(16) 



On the other hand, since Rt is by definition a real- 
positive number (\Rt\ = Rt), we can say that the dy- 
namics is given by Eq. (|13[) with E t ,Gt and Rt, defined in 
Eqs. (|I5]> and (fTgjt. 
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FIG. 2. Scheme of two initially entangled atoms, each located 
in a leaky cavity. 



III. DYNAMICS OF ENTANGLEMENT 
TRANSFER 

Having described the dynamics of the atom-cavity- 
reservoir subsystem, we now focus on the problem of two 
initially entangled atoms, each located in a leaky cavity, 
as shown in Fig. [2J This entangled state between dis- 
tant atoms have received much attention in the last years 
due to its possible applications in Quantum Communi- 
cation |11H17| . In this sense, the present study might 
be relevant. On the other hand, this coherent manipu- 
lation of atoms and cavities could also be extended into 
a trapped ions system. For example in the experimental 
setup described in reference [18|], a single trapped ion is 
coupled to a high-finesse cavity. Further improvements 
in this experiment could consider a second cavity coupled 
to a second ion. In this case, the same dynamics we are 
describing in this section can be found in this system. 
Entanglement between the electronic levels of the ions 
can be prepared using the center-of-mass mode. 

In our system, since there is no interaction between 
both atom-cavity-reservoir subsystems, the only link 
is provided by the entanglement between the atoms. 
Within this scenario, we are interested in the study of 
how the entanglement, initially shared by the two atoms, 
is transferred to other parties. To do this, we first con- 
sider an initial entangled atomic state given by 

l*o) - H.9.9)i2 + /3|ee)i 2 ) |00) 12 |00) 12 , (17) 

where cavities and reservoirs are in the vacuum state. 
Following the previous single atom-cavity-reservoir anal- 
ysis, the initial state \$o) evolves to: 



|*t) = a|s0O}i ® |.g0O} 2 + /3|^)i ® |^>2, 



(18) 



where |?/> t ) is given in Eq. (|13[) . By tracing out the degrees 
of freedom of cavities and reservoirs we are led to the 
atom-atom dynamics: 

: p 2 \E t \±\ee){ee\ + a{3\E t \ 2 (\ee)(gg\ + \gg)(ee\) 
+/3 2 \E t \ 2 (\Gt\ 2 + \Rt\ 2 )(\eg)(eg\ + \ge)(ge\) 



Pa, 



+(a 2 + p 2 (\G t \ 2 + \R t \ 2 ) \)\gg){gg\. 



(19) 



Atomic entanglement can be evaluated by using the con- 
currence [l||, leading to: 



C aia2 (i)=max{0,-2A!_ 1Q2 }, 



(20) 



where the neg ative eigenvalue of the partial transpose 
matrix [20, [2l| of p aia2 1S given by: 



A Q 



= /3\E t \ 2 ((3(l-\E t \ 2 )- 



(21) 



The entanglement flow from atoms to other parties is ob- 
tained from cavity-cavity and reservoir-reservoir reduced 
systems, respectively. The reduced cavity-cavity system 
is described by: 



PciC2 =/3 2 |G t | 4 |ll)(ll|- 
+(3 2 \G t \ 2 (\E t \ 2 



a/3|G t | 2 (|00)(ll| + |ll)(00|) 



+ (a 2 +(3 2 (\E t \ 2 + \R t \ 2 ) 



i? i | 2 )(|10)(10| 

2N2 



|oi><oil) 

)|00)(00|, (22) 



while for the reduced reservoir-reservoir system we have: 

Prir2 = /3 2 |i? t | 4 |ll)(ll| + aP\Rt\ 2 (\00)Jli\ + |11)(00|) 
+f3 2 \R t \ 2 (\E t \ 2 + \G t \ 2 ) (|10)(10| + |01)(01|) 



+ ( a 2 +/3 2 (|^| 2 + |G t | 2 ) 2 )|00)(00| 



(23) 



The corresponding entanglement in these subsystems can 
also be calculated using the concurrence. The concur- 
rences for cavities and reservoirs, respectively, are given 
by: 



C ClC2 (t) = max{0,-2Al lC2 }, 
C rir2 (t) = max{0,-2A: l, ' 2 } 



(24) 

(25) 



where the negatives eigenvalues of partial transposed ma- 
trices are given by: 



\ c l c >=p\G t \ 2 {fl{l-\G t \ 2 )-a) 
\ r 2 r >=P\R t \ 2 (p(l-\R t f)-a) 



(26) 

(27) 



We are now in a position to investigate the evolution of 
entanglement for different subsystems. Disentanglement 
dynamics depends on probability amplitudes a and /3, 
as well as the decay constant k compared to the effective 
coupling g c ff . For example, an important case is an initial 
maximally entangled state. Fig. [3] shows the entangle- 
ment evolution for atoms, cavities and reservoir subsys- 
tems, starting from this state with a = /3 for <7 ff = 5k. 
For the chosen parameters, we observe that the entan- 
glement exhibits asymptotic decay, with oscillations of 
atomic entanglement and cavities entanglement to finally 
be completely transferred into the reservoirs. These oscil- 
lations are expected because for g c g > k, there are many 
energy exchange between atomic and cavity subsystems 
before the energy is completely transferred to reservoirs. 
In addition, the entanglement between reservoirs appears 
at the same time as entanglement between atoms begins 
to decrease. This regime is reminiscent from what hap- 
pens for an initial maximally entangled state under the 
action of a dissipative environment. 

For cavities with no atoms and the initial state given by 
|f) = (c |00) + ci|ll))|00), it is well known that entan- 
glement disappears at a finite time if c\ > cq [22j . Other- 
wise, if Co > ci, the entanglement decays asymptotically. 
Moreover, in ref. [6j, finite-time disentanglement, known 
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FIG. 3. Evolution of concurrence of subsystems: a\®ai (solid 
line); c\ ® ci (dashed line) and r\ <g> r2 (dot-dashed line) for 
an initial state of the form of Eq. ()17p with a = /3 = l/\/2 
and g c g — 5k. 
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FIG. 4. Evolution of concurrence of subsystems:ai ® ai (solid 
line); c\ ® C2 (dots) and ri ® T2 (dashed line) for an initial 
state of the form of Eq. (|17|) with /3 = 1.5a and g e g — 5k. 



as, entanglement sudden death (ESD) [23|-[25j is neces- 
sarily linked to a sudden birth of entanglement between 
reservoirs, called entanglement sudden birth (ESB). This 
behavior can also occur in our system by considering un- 
balanced states with j3 > a. In Fig. HJ the evolution 
of the entanglement contained in the three subsystems 



is shown for (3 > a and g g = 5k. We realize that ei 
tanglement still experiences oscillations between caviti< 
and atoms, however, entanglement experiences sever; 
sudden deaths and entanglement sudden revivals (ESR 
whereas the reservoirs experience a sudden birth of ei 
tanglement. 

These behaviors can in principle be quantitatively ui 
derstood from Eqs. ([21]). (125)) and d^TJ). For arbitral 
effective coupling constant g c g and amplitudes a and ft, 
is not possible to calculate the times at which ESD, ES 
or ESB appears. But for the case of reservoir entangL 
ment Eq. ([27)) says that for a > f3 there will always I 
an entanglement birth at time t = 0. However, for (3 > 
the reservoirs will remain unentangled unless: 



\RtV>l--„ 



(21 



We observe this is the behavior in Fig [4] where ES 
appears for a £esb > 0. 

Some analytical calculations can be carried out in spi 
cial regimes of </ e ff as compared to k. In the strong cou- 
pling regime, when g c g 3> k we have Q w g off , so that: 



\E t \ 2 
\Gt\ 2 
\Rt\ 2 



cos {g gt)e~ 



sin 

1 - 



{g &t)e 



From Eqs. ([21)) and (|26|) we can explain the entangle- 
ment death and revivals zones for atoms and cavities in 
Fig. HJ Conditions for disentanglement in both subsys- 
tems are given by: 
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FIG. 5. Evolution of Eqs. ((32]) (solid line) and ([33]) (dashed 
line) for gr c ff = 5k. Dot-dashed line shows the value of a/ ft. 
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In the strong coupling regime, we find time windows 0.985 
with no entanglement between atoms, that is, time win- 
dows between ESD and ESR times. These time windows 
for atoms happen at time intervals that can be different q ggg 
from those for cavities, or can overlap, depending on the 
ratio a//3. At the same time the entanglement in both 
subsystems extinguish definitely because of dissipation q 975 
given by the attenuation factor in the previous inequali- 
ties. This can be observed in Fig[5] where evolutions of 
Eqs. ([52)) and (13"3l are shown compared to a/ (3. Here q g^Q 
we see that the attenuation factor exp (—Kt/2) makes it 
possible to satisfy conditions ([32]) and ([33]) . leading to 
complete disentanglement of atoms and cavities. In ad- 
dition, from Eq. (|28|) we can evaluate the time for ESB 
in the reservoir subsystem, which is given by: 
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FIG. 6. Evolution of the left-hand side of Eq. ([38]) (solid line) 
for g B a — 0.1k. Dashed line shows a//3 = 0.985. The inset 
shows the concurrence C C1C2 (£) for the cavity-cavity subsys- 
tem. 



This time is twice the found for two entangled dissipative 
cavities studied in [6]. In the present case ESB also 
happens if ft > a. 



On the other hand, in the weak coupling regime g e s <C 



k, we find that: 
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2 —Kt+^Y Kt 



(35) 



|G t | 2 ps 4 7 2 (e- 47 K * + e - Kt + 4 T Kt - 2e- Kt/2 ) (36) 



|i? t | 2 » 1 - (1 + 8j 2 )e- 4 ^ Kt + 8~f 2 e- Kt/2 , 



(37) 



where we have considered only up to second order in 
7 = ges/n- From these equations we observe that unlike 
the strong coupling regime, in the weak coupling regime 
the entanglement dynamics is not oscillatory, that is, no 
entanglement revivals can be found in atoms or cavities. 
We expect that disentanglement between atoms be fol- 
lowed by entanglement birth in reservoirs. However en- 
tanglement between cavities depends on 7 and the ratio 
a/ p. This can be understood by considering Eqs. (|26|) 
and (|36[) . such that the condition for disentanglement in 
the cavities is given by: 



1 - 4 7 2 (e- 472Kt + e - Kt+4 ~> 2,it - 2e- Kt / 2 ) > 



a 



(38) 
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The evolution of the left-hand side of the inequality is 
shown in Fig. [5] According to (|38l) , this figure shows that 
cavities are entangled depending on the ratio a/p. More 
precisely, cavities get entangled only while the curve is 
below the value of a//3. In particular, for the case of 
a/P — 0.985 (dashed line) cavities get entangled only 
when the curve is below the dashed line, as shown in 
the inset where the cavity-cavity concurrence is plotted. 
When the ratio a/P is set to a value below the min- 
imum of the left hand side in Q38p cavities never be- 
come entangled. The particular case of Fig [5] happens 
for a/P < 0.972. In such a case the entanglement ini- 
tially contained in the atomic subsystem is transferred 
directly to the reservoirs without entangling the cavities. 
This feature of entanglement dynamics can even be found 
outside the weak coupling regime as a function of 7 and 
a/ f3. We can distinguish an entangled phase and an un- 
entangled phase throughout the dynamics, which can be 
obtained from Eqs. ((24)) and p6| . Fig. [7] shows the two 
phases for the cavities entanglement dynamics as a func- 
tion of 7 and a//3. 

In the case of atom and reservoir subsystems, the non 
oscillatory behavior is shown in Fig [8] for 7 = 0.1 and 
a/P — 2/3. In such a case, the atom subsystem exhibits 
ESD, while the reservoir subsystem exhibits ESB. The 
time at which ESD and ESB happen depends on the ratio 
a//3 for a fixed 7. In the weak coupling regime, and for 
values of 7 that allows to neglect corrections in 7 2 in both 
Eqs. ((SSJ) and (|37|). we have: 
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From these Eqs. we observe that the ESB time can occur 
before, simultaneously or after the ESD. For example, in 
Fig. [SI entanglement dynamics is shown for 7 = 0.1 and 



FIG. 7. Diagram showing the entangled/unentangled phases 
of the cavity-cavity subsystem as a function of coefficients 
g c ff/n and a/p. 




FIG. 8. Evolution of concurrence of subsystems 01 ® «2, Ci (8> 
C2 and n ® r"2, using Eqs. (|35|) - (|37[ ) (solid lines) and exact 
calculations (circles). Parameter are a/j3 = 2/3 and g c g = 
O.lK. 



P = 3a. In such a figure we observe that there is a time 
window for which no entanglement is found in the three 
subsystems. Using Eqs. (|3"9"|) and (|4T)]) . we can calculate 
the size of this time window, leading to: 
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FIG. 9. Evolution of concurrence of the subsystems ai <S> a.2 
(solid line), ci Cg C2 (circles), n <S> r^ (dashed line) for /3 = 3a 
and 7 = 0.1. 



We realize from this expression that the time window 
can exist only if /3 > 2a and increases in size, as well as 
coupling strength g c g decreases. 

In this time window, where no entanglement is found 
in subsystems ai eg) 02, c\ <g> C2 and n ® r2, the question 
of where the entanglement goes in this time window be- 
comes relevant. In order to answer this question we must 
analyze the entanglement present between other bipartite 
subsystems of the overall system. 

First we consider the subsystem (01, Ci,ri) <g> 
(02,02,7*2). Since this bipartite subsystem is in a pure 
state at all times, we are able to quantify the entangle- 
ment through the square root of the tangle [261 ]. This 
entanglement is given by: 



VnaW = 2 "A 



(42) 



which corresponds to the same amount of entanglement 
initially present in the atomic subsystem. This con- 
servation of the global entanglement occurs due to the 
non- interacting character of the subsystem (01, Ci,7"i) (g) 
(a2,C2,T2) For interacting subsystems, for example, 
°i(2) ® c i(2) an d C\(2) ® r i(2)i concurrences are shown 
in Fig. QJJ From the figure, we observe that in the time 
window where no entanglement is found in subsystems 
a\ ® a2 and r\ ® r%, the interacting subsystems show 
entanglement. Indeed, in such a time window, the en- 
tanglement seems to flow only through interacting sub- 
systems. To confirm this, we need to consider different 
non- interacting subsystems such as, a\ ® C2, a\® r2 and 
c\ ®r%. In these cases the concurrences are respectively 




FIG. 10. Evolution of concurrence of the subsystems ax ® 
ci(dots), c\ eg n (circles), a\ eg) 02 (solid line) and r\ ® T2 
(dashed line) for /3 — 3a and 7 = 0.1. 
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FIG. 11. Evolution of concurrence of the subsystems 01 Cg> C2 
(dots), ai eg) ri (circles), c\ (g) r^ (squares), a\ Cg 02 (solid line) 
and n Cg r-2 (dashed line) for /3 = 3q and 7 = 0.1. 



given by: 



C aiC2 (t) = max{0, 2 [a/3E t \G t \ - y/w(Et,\G t \)] }, (43) 



C ai r 2 (t) = max{0, 2 [apEkRt - y/w(Et,Rt)] }, (44) 

C Cir2 (i) =max{0,2[a/3|G t |i? t - y/w(\G t \,Rt)]}M5) 

where w(x,y) = a i xy 2 (l — a; 2 )(l — y 2 )- The temporal 
evolution of these concurrences is shown in Fig. [IT] Ac- 
cording to this figure, in the aforementioned time win- 
dow, no contribution to entanglement comes from such 
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non-interacting subsystems. Moreover, the time win- 
dow where no entanglement is found in this subsystems 
is longer than the respective time window for the non- 
interacting subsystems a\®a,2 and r\®T2- This confirms 
that in this finite time window the entanglement can flow 
only through interacting subsystems, such as ai(2) ®C\(i) 
and Ci(2) ® r i(2): while the non- interacting subsystems 
have no entanglement. 



IV. SUMMARY 



ferent subsystems. In particular, we have found that the 
entanglement initially located in two atoms is asymptoti- 
cally mapped to the reservoirs degrees of freedom. More- 
over, although reservoirs are connected to atoms through 
the cavities, we show that for a set of initial amplitudes 
and coupling constants, there are entangled and unen- 
tangled phases in the cavities. We have also shown that 
there is a time window where no entanglement is found 
in non-interacting subsystems such us atoms, cavities or 
reservoir. In this case, we observe that the entanglement 
flows only through interacting subsystems. 



In summary, we have studied the dynamics of the en- 
tanglement transfer of two uncoupled systems each com- 
posed of a single atom inside a leaky cavity. We have 
developed a hybrid analytical approach to find the entan- 
glement dynamics without tracing out reservoir modes, 
allowing us to study the entanglement behavior in dif- 
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